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Abstract 

The article presents computer algebra procedures and routines 
applied to the study of the Dirac field on curved spacetimes. The 
main part of the procedures is devoted to the construction of Pauli 
and Dirac matrices algebra on an anholonomic orthonormal reference 
frame. Then these procedures are used to compute the Dirac equa- 
tion on curved spacetimes in a sequence of special dedicated routines. 
A comparative review of such procedures obtained for two computer 
algebra platforms (REDUCE + EXCALC and MAPLE + GRTen- 
sorll) is carried out. Applications for the calculus of Dirac equation 
on specific examples of spacetimes with or without torsion are pointed 
out. 
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1 Introduction 



In a series of recent articles [12|-|13| we have presented some routines and 



their apphcations, written in REDUCE+EXCALC computer algebra lan- 
guage for performing calculations in Dirac theory on curved spacetimes. 
Including the Dirac fields in gravitation theory requires lengthy (or cum- 
bersome) calculations which could be solved by computer algebra methods. 
Initially our main purpose was to develop a complete algebraic programming 
package for this purpose using only the REDUCE + EXCALC platform. 
This program was accomplished and the main results and applications were 
reported in our above cited articles |T^-[0. But we are aware of the fact 
that other very popular algebraic manipulations systems are on the market 
(like Mathematica or MAPLE) thus the area of people interested in alge- 
braic programming routines for Dirac equation should be much larger. In 
this perspective we developed similar programs and routines for MAPLE 
1^ platform using the package GRTensorll [|1^] (adapted for doing calcula- 



tions in General Relativity). Because there is no portability between these 
two systems we were forced to compose completely new routines, in fact fol- 
lowing the same strategy I used in REDUCE : first, the Pauli and Dirac 
matrices algebra (using only the MAPLE environment) and then the con- 
struction of the Dirac equation on curved spacetimes where the capabilities 
of GRTensorll package is used. Because the authors of GRTensorll offer also 
package versions for Mathematica we can be sure that our MAPLE routines 
could be easily adapted for Mathematica which would highly increase the 
number of users of our procedures. 

This article is organized as follows : the next section no. 2 presents a 
short review of the theory of Dirac fields on curved spacetime, pointing out 
the main notations and conventions we shall use through the article. Then 
section no. 3 is devoted to a short overview of our routines and programs 



in REDUCE+EXCALC previously described in great detail in [12]. This 
section is necessary in view of the fact that I applied the same strategy for 
constructing our programs in MAPLE as is pointed out above. Section no. 
4 contains a complete description of our programs in MAPLE+GRTensorll. 
We also included here some facts about the main differences (advantages 
and disadvantages) between the two algebraic programming platforms (RE- 
DUCE and MAPLE). Section no. 5 is devoted to the problem of including 
spacetimes with torsion in order to compute the Dirac equation using our 
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MAPLE procedures. The last section of the article includes a list of some of 
spacetimes examples we used in order to calculate the Dirac equation. Two 
of these examples are spacetimes with torsion thus it is pointed out the con- 
tribution of torsion components to the Dirac field. Several applications of 
our programs (in REDUCE or in MAPLE) are to be developed in our future 
projects : searching for inertial effects in non-inertial systems of reference 
(partially presented in |jl3| for a Schwarzschild metric without torsion) or 
quantum effects (as in |]10|) in order to provide new theoretical results for 
experimental gravity |T6|| . 

Through the article we use standard notations [|3],0, 0; Four 



dimensional spacetime indices will run from to 3 and are denoted with 
Greek letters while spatial three-dimensional indices are denoted with the 
Latin ones. The computer algebra system used was REDUCE 3.6 |^ with 
EXCALC package § or MAPLE V with GRTensorll package [T^. 



2 Pauli and Dirac matrices algebra and Dirac 
equation on curved space-time 

The main problem is to solve algebraic expressions involving the Dirac matri- 
ces , , [0 . To this end it is convenient to construct explicitly these matri- 
ces as a direct product of several pairs among the Pauli matrices ai,i = 1, 2, 3, 
and the 2x2 unit matrix. Thus all the calculations will be expressed in terms 
of the Pauli matrices and 2-dimensional Pauli spinors. Consequently the re- 
sult will be obtained in a form which is suitable for physical interpretations. 
We shall consider the Pauli matrices as abstract objects with specific multi- 
plication rules. Thus we work with operators instead of their matrices in a 
spinor representation. However, if one desires to see the result in the standard 
Dirac form with 7 matrices it will be sufficient to use a simple reconstruction 
procedure which will be presented in the next section |jl2[. 

We shall consider the Dirac formalism in the chiral form where the Dirac 
matrices are 101 : 
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The Dirac spinor : 



e H 
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(2) 



involves the Pauh spinors y?/ and (fr which transform according to the ir- 
reducible representations (1/2,0) and (0,1/2) of the group SL(2,C). In this 
representation the left and right-handed Dirac spinors are 
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(3) 



and, therefore, the Pauli spinors ipi and will be the left and the right- 
handed parts of the Dirac spinor. The SL(2,C) generators are 



[7^71 



(4) 



It is shown ^ that Tio ='H®T-L ( where 7i is the two-dimensional space of 
Pauli spinors) and, therefore the Dirac spinor can be written as: 



= il® Vl+ ^2® Vr with ^1 
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and ^2 
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while the 7-matrices and the SL(2,C) generators can be put in the form : 
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These properties allow to reduce the Dirac algebra to that of the Pauli ma- 
trices. 

In order to introduce Dirac equation on curved space-time we have always 
used an anholonomic orthonormal frame because at any point of spacetime 
we need an orthonormal reference frame in order to describe the spinor field 
as it is already pointed before [l^). The Dirac equation in a general reference 
of frame, defined by an anholonomic tetrad field is |^ : 

ikj^D^^H = mc^' (8) 



where the covariant Dirac derivative is 



and where S'^'^ are the SL(2,C) generators 
components of the connection. 



and r 



upfj. 



(9) 

are the anholonomic 
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3 Review of the REDUCE+EXCALC rou- 
tines for calculating the Dirac equation 

We shall describe here those part of the program realizing the Pauli and 
Dirac matrix algebras |12|,[0. In the first lines of this sequence we introduce 
the operators and the non-commuting operators being useful throughout the 
entire program. The Pauli matrices are represented using the operator p with 
one argument. The Dirac matrices are denoted by gam of one argument (an 
operator if we use only REDUCE, or for EXCALC package it will be a 0-form 
with one index) while the operator dirac stands for the Dirac equation. The 
SL(2,C) generators are denoted by the 0-form s(a,b).The basic algebraic 
operation, the commutator (com) and anticommutator (acom) are then 
defined here only for commuting (or anticommuting) only simple objects 
("kernels"). For commuting more complex expressions, (in order to introduce 
some necessary commutation relations) a more complex operator is necessary 
to introduce. Other objects, having a more or less local utilization in the 
program will be introduced with declarations and statements at their specific 
appearance. 

The main part of the program is the Pauli subroutine : 
LET p(0)=l; 

LET p(2)*p(l)=-p(l)*p(2) 



LET p(l)*p(2)=i*p(3) ; 
LET p(3)*p(l)=-p(l)*p( 
LET p(l)*p(3)=-i*p(2); 
LET p(3)*p(2)=-p(2)*p( 
LET p(2)*p(3)=i*p(l) ; 
LET p(l)**2=l 
LET p(2)**2=l 
LET p(3)**2=l 

The Pauli matrices, cTj appear as p(i) while the 2x2 unity matrix is p(0)=l. 
The properties of the Pauli matrices are given by the above sequence of 10 
lines. The direct product denoted by the pd operator has properties, 
introduced as : 

for all a,b,c,u let pd(a,b) *pd(c ,u)=pd(a*c ,b*u) ; 
for all a,b let pd(a,b) **2=pd(a**2 ,b**2) ; 
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for all a,b let pd(-a,b)=-pd(a,b) ; 
for all a,b let pd(a,-b)=-pd(a,b) ; 
for all a,b let pd(i*a,b)=i*pd(a,b) ; 
for all a,b let pd(a, i*b)=i*pd(a,b) ; 
for all a let pd(0,a)=0; 
for all a let pd(a,0)=0; 
let pd(l,l)=l; 

Some difficulties arise from the bilinearity of the direct product which requires 
to identify all the scalars involved in the current calculations. This can be 
done only by using complicated procedures or special assignments. For this 
reason we shall use a special definition of the direct product (pd) which 
gives up the general bilinearity property. The operator pd will depend on 
two Pauli matrices or on the Pauli matrices with factors —1 or ±i. It is able 
to recognize only these numbers but this is enough since the multiplication 
of the Pauli matrices has just the structure constants ±1 and ±i (we have 
ataj = 6ij + ieijkCTk). 

Thus by introducing the multiplication rules of the direct product it will 
be sufficient to give some instructions (see above) which represent the bi- 
linearity, defined only for the scalars —1 and i. The next two instructions 
represent the definition of "0" in the direct product space, while the last one 
from the above sequence introduces the 4x4 unit matrix. 

Thus the 7-matrices can be defined now with the help of our direct prod- 
uct; we added also here the definition of the SL(2,C) generators from eq. (^) 

gamd) :=i*pd(p(2) ,p(l)) ; 
gam(2) :=i*pd(p(2) ,p(2)) ; 

gam(3) :=i*pd(p(2) ,p(3) ) ; % Remember that gam's are 

gain(0):= pd(p(l),l); 7„ 0-forms !!! 

gamS := -pd(p(3) , 1) ; % instead of "gam(5) " 

s(a,b):= i*com(gam(a) ,gam(b) )/4; 

All the above program lines we have presented here can be used for the 
Dirac theory on the Minkowski space-time in an inertial system of reference. 
Here we shall point at first the main differences which appear when one 
wants to run our procedures on some curved space-times or in a non-inertial 
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reference of frame. Some of these minor differences are already integrated in 
the hnes presented in the previous section. 

First of all we have to add, at the beginning of the program some EX- 
CALC lines containing the metric statement. We have always used an an- 
holonomic orthonormal frame because at any point of spacetime we need an 
orthonormal reference frame in order to describe the spinor field as it has 
been already pointed before [ll2| ). 

After the above metric statement we must add in the program all the 
procedures described in the last section. Then we can introduce some lines 
calculating the Dirac equation in this context. As a result we have used the 
next sequence of EXCALC lines : 

pform {der(j) ,psi}=0; fdomain psi=psi (x,y,z,t) ; 
der(-j):= ee(-j)_|d psi + (i/2)*s(b,h)*cris(-j ,-b,-h) ; 

operator derpO , derpl , derp2 , derp3 ; 
noncom derpO, derpl ,derp2,derp3; 

let (§(psi ,t)=derpO ; 
let (§(psi ,x)=derpl ; 
let (§(psi ,y)=derp2 ; 
let (§(psi ,z)=derp3; 

dirac := i*has*gam(j)*der(-j)-m*c; 

hajii:= -(gam(0)*(l/(ee(-0)_|d t) ) *dirac-i*has*derpO) ; 

In defining the Dirac derivative der we have introduced also an formal Dirac 
spinor (psi) being a 0-form and depending on the variables imposed by the 
symmetry of the problem. It is just an intermediate step (in fact a "trick") in 
order to obtain the partial derivative components as operators, because after 
calculating the components of the covariant derivative (der(-j)- see above) 
we have to replace the partial derivatives of psi with four non-commuting 
operators derpO, derpl ... derp3. The Dirac operator is thus defined 
as dirac: = [ikj^^D^ — mc)\E' and finally the Dirac Hamiltonian (ham) is 
obtained from the canonical form of the Dirac equation : 

ih— = Hip 
at 
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which we shall use later, in the study of the nonrelativistic approximation of 
the Dirac equation in non-inertial reference frames. 

The results we have obtained after processing the program lines presented 
until now contain only the Pauli matrices and direct products of Pauli ma- 
trices. If one wants to have the final result in terms of the 7-matrices and 
SL(2,C) generators (and not in terms of direct products of Pauli matrices) 
the procedure rec should be used : 

operator gama,gen; 
noncom gama,gen; 
PROCEDURE rec(a) ; 
begin; 

wsl :=sub(pd(p(l) , l)=gama(0) ,a) ; 
wsl :=sub(pd(p(2) , l)=-i*gama(0) *gama(5) ,wsl) ; 
wsl :=sub(pd(p(3) , l)=-gama(5) ,wsl) ; 
wsl:=sub(pd(l,p(l))=2*gen(2,3) ,wsl) ; 
wsl:=sub(pd(l,p(2))=2*gen(3,l) ,wsl) ; 
wsl:=sub(pd(l,p(3))=2*gen(l,2) ,wsl) ; 
for k:=l:3 do 

«wsl :=sub(pd(p(2) ,p(k) )=-i*gama(k) ,wsl) ; 

wsl :=sub(pd(p(l) ,p(k))=gama(k)*gama(5) ,wsl) ; 
wsl:=sub(pd(p(3) ,p(k))=2*i*gen(0,k) ,wsl)»; 
return wsl; 
end; 

This is an operator depending on an expression involving matrices (a) which 
reconstructs the 7-matrices and the SL(2,C) generators from the direct prod- 
ucts of Pauli matrices according to eqs. (6) and (7). 

We must point out that the new introduced operators gama and gen 
do not represent a complete algebra. They are introduced in order to have 
a result in a comprehensible form. Thus, in this form the result cannot be 
used in further computations. Only the results obtained before processing 
the rec procedure could be used, in order to benefit of the complete Pauli 
and Dirac matrices algebra. 
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4 MAPLE+GRTensorll procedures for cal- 
culating the Dirac equation 

Here we shall present, in details, our procedures in MAPLE+GRTensorll 
for calculating the Dirac equations, pointing out the main differences be- 
tween MAPLE and REDUCE programming in obtaining the same results. 
The first major problem appears in MAPLE when one try to introduce the 
Pauli and Dirac matrices algebra. In MAPLE this will be a difficult task 
because the ordinary product (assigned in MAPLE with "*" ) of operators is 
automatically commutative, associative, linear, etc. like an ordinary scalar 
product - in REDUCE these properties are active only if the operators are 
declared previously as having such properties. Thus we have to define two 
special product operators : for Pauli matrices cTq,, a = 0, 1, 2, 3 (assigned in 
our procedures with pr) and for the direct product of Pauli matrices (as- 
signed here also with the operator pd(",")) which is assigned with "&;p". 
As a consequence we have to introduce long lists with their properties as, for 
example : 

> def ineCsigma, sigma(0)=l) ; 

> def ine(pr,pr(l , 1)=1 ,pr(l ,sigma(l) )=sigma(l) ,pr (1, sigma(2))= 
sigma(2) ,pr(l,sigma(3))=sigma(3) ,pr (sigma(l) ,l)=sigma(l) , 
pr(sigma(2) , l)=sigma(2) ,pr(sigma(3) , l)=sigma(3) , . . . 
pr(sigma(l) ,sigma(l))=l,pr(sigina(2) , sigma(2) )=1 , . . . 

> def ine(pd,pd(0,a: : algebraic)=0,pd(a: : algebraic ,0)=0,pd(l , 1)=1 , 
pd(I*a: :algebraic,b: : algebraic)=I*pd(a,b) , 

pd(-I*a: : algebraic, b: : algebraic)=-I*pd(a,b) , 
pd(a: : algebraic, I *b: : algebraic)=I*pd(a,b) , 
pd(a: : algebraic, -I*b: : algebraic)=-I*pd(a,b) , 
pd(-a: : algebraic ,b : : algebraic) =-pd (a, b) , 
pd(a: : algebraic, -b: : algebraic) =-pd (a, b) ) ; 

> def ineC&p' , '&p' (-a: : algebraic, -b: : algebraic) = '&p' (a,b) , . . . 

Of course the reader may ask why is not much simpler to declare, as an 
example, the &:p as being linear (or multilinear) ? Because in this case the 
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operator does not act properly, the linearity property picking out from the 
operator all the terms, being or not Pauli matrices or direct products pd of 
Pauli matrices. Thus it is necessary to forget the linearity and to introduce, 
as separate properties all the possible situations to appear in the calculus. 
As a result the program becomes very large with a corresponding waste of 
RAM memory and speed of running. This will be the main disadvantage of 
MAPLE version of our program in comparison with the short (and, why not, 
elegant) REDUCE procedures. Of course, in a more compact version of our 
programs, we defined MAPLE routines with these operators, and the user 
needs only to load at the beginning of MAPLE session these routines, but 
there is no significant economy of memory and running time. 

The next step is to define Dirac 7-matrices and a special commutator 
(with &;p) : 

> define (gam, 

gain(l)=I*pd(sigma(2) ,sigma(l)) , 
gam(2)=I*pd(sigma(2) ,sigma(2)) , 
gam(3)=I*pd(sigma(2) ,sigma(3)) , 
gam(0)=pd(sigma(l) ,1) ,gam(5)=-pd(sigma(3) ,1)) ; 

> def ine(comu,comu(a: :algebraic,b: :algebraic)=a &p b - b &p a); 

The next program-lines are in GRTensorll environment, thus is necessary 
to load first the GRTensorll package and then the corresponding metric (with 
qload(...) command |15|). It follows then : 

> grdef('SS{ "a "b }'); 

> grcalc(SS(up,up) ) : 

> (I/4)*comu(gaiii(0) ,gain(0)) ; 

> (I/4)*comu(gaiii(l) ,gain(0)) ; 



> (I/4)*comu(gaiii(3) ,gain(3)) ; 

> grdisplay(SS(up,up)) ; 

> grcalc(Chr(dn,dn,dn)) ;grdisplay(Chr(dn,dn,dn)) ; 

> grcalc (Chr (bdn , bdn , bdn) ) ; grdisplay (Chr (bdn , bdn , bdn) ) ; 
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These are a sequence of commands in GRTensorll for defining the SL(2,C) 
generators S^^ (as the tensor SS{ 'a 'b}) using formula and for the cal- 
culus of Christoffel symbols in an orthonormal frame base (Chr(bdn,bdn)). 
Here it becomes obvious one of the main advantages of MAPLE+GRTensorll 
platform, namely, the possibility of computing of the tensor components both 
in a general reference frame or in an anholonomic orthonormal frame which 
is vital for our purpose of construction the Dirac equation. 

Next we have to define, as two vectors the Dirac-7 matrices (assigned as 
the contravariant vector ga{ 'a } and the derivatives of the wave function 
ip (assigned as the covariant vector Psid{ a } in order to use the facilities of 
GRTensorll for manipulating with indices : 

> grdef('ga{ "a }: = [gain(0) ,gaiii(l) ,gam(2) ,gain(3)] ') ; 

> grdisplay(ga(up)) ; 

> grdef('Psid{ a }:=[diff (psi(x,t) ,t) ,diff (psi(x,t) ,x) , 

diff (psi(x,t) ,y) ,diff (psi(x,t) ,z)] ') ; 

> grcalc(Psid(dn) ) ;grdisplay(Psid(dn) ) ; 

> grcalc(Psid(bdn) ) ;grdisplay(Psid(bdn) ) ; 

The next step is to define the term i^SP^Vypi^, from equation (|^) : 

> grdef('de{ i }:=(I/2)*SS{ "a "b }*Chr{ (i) (a) (b) }'); 

> grcalc(de(dn)) ;grdisplay(de(dn) ) ; 

The components of de{ i } are polynomials containing direct products 
pd(...) of Pauli matrices and the product between Y ^"^^ ^S^^T^p^ from 
equation (||) is, in fact, the special product &ip. Thus we obtain the term 
■y'^S^'^T^pp (denoted below with the operator dd) by a special MAPLE se- 
quence which in fact split the components of de{ i } in monomial terms and 
then execute the corresponding 8zp product, finally reconstructing the dd 
operator : 

aO : =expand(gr component (de (dn) , [t] ) ) ; aOO : =0 ; 

> uO : =whattype (aO) ; uO ; nops (aO) ; 

> if uO='+' then for i from 1 to nops(aO) do 

aOO : =aOO+I*h*grcomponent (ga(up) , [t] ) &p op(i,aO) od else 
aOO : =I*h*grcomponent (ga(up) , [t] ) &p aO fi; aOO; 
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> al : =expand (grcomponent (de (dn) , [x] ) ) ; al 1 : =0 ; 

> ul :=whattype(al) ;ul; 

> nops(al) ; 

> if ul='+' then for i from 1 to nops(al) do 

all : =all+I*h*grcomponent (ga(up) , [x] ) &p op(i,al) od else 
all : =I*h*grcomponent (ga(up) , [x] ) &p al fi;all; 



> grdef Cdd') ; grcalc(dd); 

> a00+all+a22+a33; 

> grdisplay (dd) ; 

Finally the Dirac equation is obtained as : 

> grdef (' dirac :=I*h*ga{ '"l }*Psid{ (1) } + dd*psi(x,t) - m*c*psi (x,t) ' ) ; 

> grcalc (dirac) ; 

> grdisplay (dirac) ; 

In order to obtain the Dirac equation in a more comprehensible form we have 
the next sequence of MAPLE commands (similar to the reconstruction rec 
procedure from the REDUCE program) : 

> def ine( 'gen' ) ; 

> define ( 'gama' ) ; 

> grmap (dirac , subs ,pd(sigma(l) , l)=gama(0) ,pd(sigma(2) , 1)=-I*gama(0)* 
gama(5) ,pd(sigma(3) ,l)=-gama(5) ,pd(l,sigma(l))=2*gen(2,3) ,pd(l, 
sigma(2))=2*gen(3,l) ,pd(l,sigma(3))=2*gen(l,2) ,pd(sigma(2) , sigma(l) )= 
-I*gama(l) ,pd(sigma(2) , sigma(2) )=-I*gama(2) ,pd(sigma(2) ,sigina(3))= 
-I*gama(3) ,pd(sigma(l) ,sigma(l))=gama(l)*gama(5) ,pd(sigma(l) ,sigma(2))= 
gama(2) *gama(5) ,pd(sigma(l) , sigma(3) )=gama(3) *gajna(5) ,pd(sigma(3) , 
sigma(l))=2*I*gen(0,l) ,pd(sigma(3) ,sigma(2))=2*I*gen(0,2) ,pd(sigma(3) , 
sigma(3))=2*I*gen(0,3) , 'x') ; 

> grdisplay (dirac) ; 

where, of course, as in the REDUCE version, the operators gen and gama 
do not represent a complete algebra. 
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5 Dirac equation on spacetimes with torsion 



We shall present here the way we adapted our MAPLE+GRTensorll pro- 
grams in order to calculate the Dirac equation on space-times with torsion. 

The geometrical frame for General Relativity is a Riemannian space- 
time but one very promising generalization is the Riemann-Cartan geometry 
which (i) is the most natural generalization of a Riemannian geometry by 
allowing a non-symmetric metric-compatible connection, (ii) treats spin on 
the same level as mass as it is indicated by the group theoretical analysis 
of the Poincare group, and (iii) arises in most gauge theoretical approaches 
to General Relativity, as e.g. in the Poincare-gauge theory or supergravity 
However, till now there is no experimental evidence for torsion. On 
the other hand, from the lack of effects which may be due to torsion one can 



calculate estimates on the maximal strength of the torsion fields [|^. In this 
aspect we think that it is possible, using computer algebra facilities to ap- 
proach new theoretical aspects on matter fields (for example the Dirac field) 
behavior on spacetimes with torsion in order to point out new gravitational 
effects and experiments at microscopic level. 

A metric compatible connection in a Riemann-Cartan theory is related 
to the torsion components by (see |]^ - eq. (1.18)) 

ra/37 = f — - [(Ca/37 — Cp^a + C-yajs) — {Taf3-y — Tp^a + ^70/?)] (10) 

where Tap-y are the components of the riemannian connection, Cafi-y is the 
object of anholonomicity and T^/s-y are the components of the torsion. The 
idea is to replace the connection components from the covariant derivative 
appearing in eqs. (|§-|^) with the above ones, of course after calculating them 
in an orthonormal anholonomic reference frame, suitable for calculation of 
the Dirac equation. Thus the sequence for calculating the de{dn} operator 
(see above) should be replaced with 

> grdef('de{ i }:=(I/2)*SS{ "a "b }*(CHR{ (i) (a) (b) }+ 

(l/2)*(tor{ (i) (a) (b) } -tor{ (a) (b) (i) } 
+ tor{ (b) (i) (a) })))'); 

> grcalc(de(dn) ) ; grdisplay(de(dn)) ; 

where the new connection components CHR{a, b, c} are now defined by the 
sequence 
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> grdef('ee-[ a "b }:= wl{ "b }*kdelta{ a $x } 

+ w2{ "b }*kdelta{ a $y } + 
w3-[ "b }*kdelta{ a $z } + w4{ "b }*kdelta-[ a $t }'); 

> grcalc(ee(dn,up)) ; 

> grdisplay (ee (dn,up) ) ; 

> grdef('CC-[ a b c }:=2*ee{ a ~i }*ee-[ b ~ j }*ee{ c [j ,i] }'); 

> grcalc(CC(dn,dn,dn)) ; grdisplay (CC (dn, dn, dn) ) ; 

> grdef('CHR{ (a) (b) (c) } := Chr{ (a) (b) (c) } - 

(1/2)*(CC{ a b c } - 
CC{ b c a } + CC{c a b }) ') ; 

> grcalc(CHR(bdn,bdn,bdn)) ; grdisplay (CHR(bdn,bdn,bdn) ) ; 

and the rest of the routines are unchanged. The only problem remains now 
to introduce in an adequate way the components of the torsion tensor. I used 
the suggestion from ^ , pointing that we can assume that the 2- form T° asso- 
ciated to the torsion tensor should have the same pattern as the dd'^^s where 
6"' is the orthonormal coframe, who's components are denoted in GRTen- 
sorll with wl{dn} ...w4{dn}. Thus this operation is possible only after we 
introduced the metric (with qload command). Then calculating the deriva- 
tives of the orthonormal frame vector basis components we can introduce the 
torsion components by inspecting carefully these derivatives. Here there is 
an example of how this becomes possible in MAPLE+GRTensorll using one 
of the metric examples presented in the next section: 

> grcalc (wl (dn,pdn) ) ;grcalc(w2(dn,pdn) ) ; 

> grcalc (w3(dn,pdn) ) ; grcalc (w4(dn,pdn) ) ; 

> grdisplay (wl(dn,pdn) ) ; grdisplay (w2 (dn, pdn) ) ; 

> grdisplay (w3(dn, pdn) ) ; grdisplay (w4(dn, pdn) ) ; 

> grcalc (wl (bdn , pbdn) ) ; grcalc (w2 (bdn , pbdn) ) ; 

> grcalc (w3 (bdn , pbdn) ) ; grcalc (w4 (bdn , pbdn) ) ; 

> grdisplay (wl(bdn, pbdn) ) ; grdisplay (w2 (bdn, pbdn) ) ; 

> grdi splay (w3 (bdn, pbdn) ) ; grdisplay (w4 (bdn, pbdn) ) ; 

> grdef ('tor{ "a b c }:=wl{ b ,c }*kdelta-[ "a $x } + 

w2{ b ,c }*kdelta{ "a $y } + 
w3{ b ,c }*kdelta{ "a $z } + 
w4{ b ,c }*kdelta{ "a $t }'); 

> grcalc (tor (up, dn, dn) ) ; grdi splay (tor (up, dn,dn)) ; 

> grmap(tor(up,dn,dn) ,subs,f (x)=v4(x) ,h(x)=v3(x) ,g(x)=v2(x) , 'x') ; 
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> grdi splay (tor (up, dn, dn) ) ; 

> grcalc(tor(bup,bdn,bdn) ) ; 

> grdisplay (tor (bup , bdn , bdn) ) ; 

The reader can observe that we first assigned the components of the torsion 
tensor (here denoted with tor{up,dn,dn}) then after displaying his com- 
ponents we can decide to substitute new functions describing the torsion 
instead of the functions describing the metric. Of course finally we calcu- 
late the components of the torsion in an orthonormal anholonomic reference 
frame (tor{bup,bdn,bd}). 



6 Some specific results 

This section is devoted to a list of some recent results we obtained by running 
our procedures in MAPLE+GRTensorll, already described in the previous 
sections. First we tested our programs by re-obtaining the form of the Dirac 
equations in several spacetime metrics, obtained with REDUCE+EXCALC 



procedures and reported in our previous articles [|T3]-[]T3|- The concordance 
of these results with previous ones was a good sign for us to proceed with 
more complicated and new examples, including ones with torsion. Here we 
shall present some of these late examples. 



1. Conformally static metric with $ and S constant where the line element 
is : 

^^2 ^ g2$<+2E^^^)2^^2 ^ g2<I.t+2S^2^^2 ^ g2*t+2E^2 5,^(^)2^^2 _ g2$t+2E^(^)2^^2 

Thus the Dirac equation becomes : 



ihe 



-<S>t-T, 



1 



1 d 



^a{r) dt 2a{r)b{r) 



^ h'ir) - 



a[r)r ^ 
1 



1 d 



+ 



d 



1 1 n/3 d 



r dr r sin{6) 06 
ip — mcip = 



where b'{r) is the derivative db{r)/dr. 
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2. Tauh-NUT spacetime having the hne element as : 

ds^ = -AfU{t)dy^ - 8/2f/(t) cos{e)dyd(j) - {f + l^)de'^ 
+ {-4fU{t) cos{ef - {e + sin{ef)d4? + 



U{t) 



df 



the coordinates being {y,9,(j),t). We obtained the Dirac equation as 



ih 



t^ + P 



t^ + f 



-7 



U\t) + JU{t) 1 + 



dt 



cotg{e)VWTP-f 



m 

mcipit) = 



3. Gddel spacetime, having the hne element as : 

ds^ = -a^dx^ + ^a^e^^dy^ + 2a^e''cdydt - a^dz^ + a^c^dt^ 

in coordinates {x,y,z,t). Here the Dirac equation is simply : 
d \ n: o ( ^ d 1 d\ ,(9 Id 



ih— 
a 



^ I 2 dxj ^ \ dy cdtj ^ dz cdt 



^p{x,y,z,t) 
—mcipix, y, z,t) = 
4. McCrea static spacetime with torsion having the line element as 

in (x, y, z, t) coordinates. If the coordinate lines of y are closed with < 
y < 271 and — oo <2;<oo, 0<a;<oo, the spacetime is cylindrically with 
y as the angular, x the cylindrical radial and z the longitudinal coordinate. 
If —oo,x,y,z < oo the symmetry is pseudo-planar. In McCrea considers 
the simplest solution of Einstein equations with cosmological constant as 

f = h = h = qx/3 (11) 

and the cosmological constant turns to be We shall first consider the 

general case specializing the results at the final step of the program to the 
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above particular solution. Running our MAPLE+GRTensorll procedures, at 
first we obtain the torsion tensor component as : 



the rest of the components being zero. This time we have obtained the Dirac 
equation, depending also on the components of the torsion tensor as : 

^/^i^i Lv2(x) dv2{x) ^ ^^s(^) dv3{x) ^ ^^i(^) dv4:{x) _ dfjx) _ dg{x) 
2 \ dx dx dx dx dx 

Qfiix) \ 

^ ' + 2— ^[x) - mc\l)[x) = 



dx dx ^ 

Of course when we specialize to the particular solution proposed by McCrea 
in [PI we have to assign the form of metric functions as in (|lT]) and we can 
then take the torsion functions as 

v2 = v3 = vA = v{x) 

and the Dirac equation becomes 

5 Schwarzschild metric with torsion. This example is interesting in the view 
of recent investigations on the contribution of torsion in gravity experiments 



using atomic interferometry \U]\. Here we used the Schwarzschild metric 



having the line element written as 

Here we prefer to specialize the form of A(r) and z/(r) functions as 

2m 1 



z/(r) = 1 



A(r) 



(we have G = c = 1) later, after obtaining the form of the Dirac equation in 
term of A(r) and v{r) functions. 
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Using the same "trick" as in the previous example, we choose the com- 
ponents of the torsion tensor as 

= ^^eV^/^C^) ; Tl = 1 ; = sin{e) 



Tl = ^^eV^/^(0 ; = r cos{e) 



the rest of the components being zero. Running away our procedures we 
obtain the Dirac equation containing terms with torsion tensor components 
as : 



ih, 



1 -A(.)/2^1 (^[1, 



7 



,/2(r)/2 _ 



dr 



l^\cos{d) - cotg{e)) 
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